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Abstract 

We consider a partially ionized gas at thermal equilibrium, in the Saha regime. 
The system is described in terms of a quantum plasma of nuclei and electrons. In this 
framework, the Coulomb interaction is the source of a large variety of phenomena oc- 
curing at different scales: recombination, screening, diffraction, etc. In this paper, we 
derive a cluster expansion adequate for a coherent treatment of those phenomena. The 
expansion is obtained by combining the path integral representation of the quantum 
gas with familiar Mayer diagrammatics. In this formalism, graphs have a clear phys- 
ical interpretation: vertices are associated with recombined chemical species, while 
bonds describe their mutual interactions. The diagrammatical rules account exactly 
for all effects in the medium. Applications to thermodynamics, van der Waals forces 
and dielectric versus conductive behaviour will be presented in forthcoming papers. 

Key words: Ionized gases; quantum plasmas; cluster expansions; recombination; screen- 
ing 

1 Introduction 

This paper is devoted to the construction of a cluster expansion adapted to the study of 
partially ionized gases. In the so-called physical picture of such gases, one considers an 
assembly of non relativistic point nuclei and electrons in thermal equilibrium with appro- 
priate statistics and interacting with the Coulomb potential. That fundamental approach 
is to be contrasted with the chemical picture which introduces preformed chemical species 
as elementary constituents. At this basic level of description all characteristic phenom- 
ena of charged systems such as atomic or molecular recombination, ionization, ionic and 
dielectric screening originate from the sole Coulomb interaction between charges. It is 
notoriously difficult to treat all of them from first principles in a fully coherent way. In 
particular one has to deal simultaneously and without intermediate modelization with the 
short distance part of Coulomb potential which gives rise to chemical binding (at the scale 
of Bohr radius), and with its long range part responsible for collective screening effects (at 
the scale of the classical Debye length). 
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An important point should be stressed straight away. The very concept of atom and 
molecule can only make sense in a low density or moderately dense medium and at tem- 
peratures sufficiently low to be away from full ionization. This situation is captured in a 
scaling limit called the atomic (or molecular) limit, which is more precisely described in 
Section 5 in the case of the Hydrogen plasma. In that limit, one lets temperature and 
density tend to zero in a coupled way. Low temperature favors binding over ionization, 
whereas low density, by increasing the available phase space, favors dissociation. The 
rate at which density is reduced as the temperature tends to zero determines an energy 
entropy balance that selects the formation of some chemical species. In the strict limit 
one recovers the thermodynamical laws of mixtures of ideal substances. Thus the proper 
physics of partially ionized gases lies in the vicinity of the atomic limit: this defines the 
Saha regime for which it is of interest to provide a suitable expansion. 

In the literature, there exists various attempts, within the physical picture, to de- 
rive suitable expansions at low temperature and low density. First, the effective-potential 
method has been applied to quantum plasmas by Ebeling 2 J. That first-principles ap- 
proach amounts to introduce an equivalent classical system of point particles with many- 
body potentials which incorporate all quantum effects. Low density expansions at fixed 
non-zero temperature have been successfully computed while plausible expres- 

sions for basic ingredients of the chemical picture (like effective fugacities of recombined 
entities) have also been proposed within that framework OEl- Unfortunately, an explicit 
control of the contribution arising from any n-body potential (n > 3) cannot be achieved 
by using simple diagrammatical techniques (in most previous works, only two-body po- 
tentials are retained). Therefore, all effects related to recombination into complex entities 
of more than two particles, cannot be properly taken into account. Another approach jS] 
relies on a heuristic introduction of quantum and screening effects in classical versions of 
usual fugacities expansions. Roughly speaking, Boltzmann factors at short distances are 
replaced by their quantum counterparts, while Coulomb interactions at large distances are 
exponentially screened over classical Debye screening length. If that method may provide 
rather good quantitative results at moderate temperatures and densities [Jj, explicit ex- 
pressions for contributions related to any phenomenon are not available. Eventually, we 
emphasize that many-body perturbation theory [SJ |§1 is not suited to the description of 
low temperature and low density regimes. Indeed, perturbative expansions with respect 
to the coupling constant (the electric charge) cannot account for recombination where 
interactions play a crucial non-perturbative role. Consequently, that formalism is more 
appropriate for dealing with fully ionized regimes, either at low temperature and high 
density ^01 , or at high temperature and low density 11 2 j . 

In order to circumvent the drawbacks of previous methods, we introduce the path 
integral representation of quantum Coulomb gases. As shown in this paper, that formalism 
allows us to build an expansion which fulfills all requirements for a proper description of 
Saha regime. In particular, collective screening effects, as well as quantum mechanics of 
arbitrary particle clusters, are simultaneously dealt with. For the sake of presentation 
simplicity, derivations are carried out for the Hydrogen plasma made with protons and 
electrons (see Section 2). Extensions to other quantum gases are briefly mentionned in 
Section 5. 

The functional integral representation of quantum Coulomb gas by means of Feynman- 
Kac formula 13 is displayed in Section 3.1. Charges of the same kind and same statistics 
are collected according to the cycle decomposition of the permutation group into closed 
Brownian paths called charged loops |14| . In terms of those loops, statistical averages 
are performed according to the same rules as in classical statistical mechanics. As a 
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consequence, the powerful methods of Mayer graphs are available in the space of loops 
|15j : activities including the self-energy of a loop are associated to vertices, and bonds are 
usual Mayer factors of interaction between different loops. Because of the long range of 
the Coulomb potential, those bonds are not integrable, and one has first to resum classes 
of diverging graphs in order to obtain an effective potential with shorter range. Since 
details have already been worked out in previous works |161 115j we merely summarize the 
procedure in section 3.2 for the sake of completeness. The procedure follows the familiar 
chain graph summation as in the classical case ^1 The subsequent effective 

potential is the quantum analogue of the Debye potential. It incorporates specific features 
of quantum mechanics, namely quantum statistics, short distance Gibbs weight leading 
to binding, as well as algebraically decaying tails due to quantum fluctuations. That 
effective potential is extensively studied in |2Qj: with its help one can define a resummed 
diagrammatics of prototype graphs that are free from Coulomb divergencies. The rules 
are the same as those of standard Abe-Meeron diagrammatics. This fugacity expansion 
in terms of loop prototype graphs is suited to derive, at fixed non-zero temperature, low 
density virial equation of state in the ionized phase |21[ I22j . The small parameter is the 
density p e of unbound electrons so that contributions of bound states (e.g. hydrogen 
atoms) only occur as small corrections at order f? e . 

The expansion (|3.25|) in terms of loop prototype graphs is not well adapted to a direct 
analysis of Saha regime, where recombined atoms or molecules appear with comparable 
densities at low temperature. The reason for this inadequacy is that Gibbs factors arising 
in a given prototype graph do not provide the complete set of pair interactions that would 
correspond to an atomic or molecular Hamiltonian. Hence each prototype graph is not in 
direct relation with the quantities that are naturally involved in the definitions of energies 
and densities of recombined entities. Then, the aim is to systematically reorganize the 
previous loop expansion into an expansion where clusters of charges are fully interacting 
and are in one to one correspondance with the possible chemical binding processes. This 
is the content of Section 4, where we derive the main result of the paper, namely new form 
(|4.22j) of the fugacity expansion that we call the screened cluster expansion. 

The derivation is accomplished in three steps. The first one (Section 4.1) consists 
in collecting prototype graphs into classes such that the sum of graphs in a class builds 
weights of fully interacting clusters of loops. This provides Gibbs factors where all pairwise 
interactions between loops are present. The building rules of those coefficients are the same 
as those of usual Mayer cofficients in standard fugacity expansions. In fact, had one to deal 
with a gas of particles interacting by means of short range forces, those coefficients would 
simply be identical to standard Mayer coefficients written within Feynman-Kac functional 
integral formalism. The difference here is that loops interact with the effective potential 
requested by the long range of Coulomb interaction. As that effective potential itself results 
from chain summations, one must be careful to avoid double counting of graphs. It turns 
out that avoiding double counting induces residual interactions between clusters through 
the effective potential. If now Mayer coefficients of fully interacting clusters are interpreted 
as dressed activities associated with vertices and residual interactions are associated with 
bonds, it is nice that the resulting diagrammatics is still of Mayer type in terms of those 
weights and bonds: this leads to formally exact representation q4.5|) of the loop density. 
At any non-zero temperature, each graph of (|4,5j) is finite since bonds and vertices are 
entirely constructed with the help of the effective potential. 

The necessity of the second step performed in Section 4.2 is the following. As explained 
before, the physical range of thermodynamical parameters relevant for the description of 
partially ionized gases lies in the neighborhood of the atomic limit, namely for vanishing 
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densities and temperatures. If one lets the density tend to zero in weights of (|4,5j) . the 
effective potential reduces to the bare Coulomb potential so that one faces again diverging 
spatial integrals. Hence one needs to introduce a further step if one wants cluster integrals 
to have a well defined atomic limit. In rough terms, it consists simply in a regularization 
procedure of divergent integrals by substraction of non-integrable part of their integrands. 
This leads to the introduction of truncated weights in such a way that new weights remain 
integrable when screening is removed. The definition of truncated cluster weights adopted 
here has a remarkable property: it preserves the topological rules of Mayer diagrammatics. 
This leads to expansion (|4.18j) for the loop density. The vertices are truncated Mayer 
coefficients, and a bond between two clusters links every loop in one cluster to every loop 
in the other cluster via the effective potential (or its second and third power). 

In a third step (Section 4.3), we introduce particle clusters. If in each graph of (|4.18|) . 
we detail phase-space integrations over loop internal degrees of freedom, loop clusters 
give rise to particle clusters. The final step amounts to collect all contributions with the 
same topological structure in terms of particle clusters. Then, we obtain the so-called 
screened cluster expansion (|4.22|) for the particle density. In that diagrammatic series, 
the graphs G (see Fig. 12) are identical to familiar Mayer graphs with particle clusters in 
place of one-particle points. The internal state of a cluster is determined by a distribution 
of its particles into a given set of loops. Its statistical weight reduces to the truncated 
Mayer coefficient for that loop configuration. Two clusters are connected by at most one 
bond, which can be either the total screened interaction between the respective two sets 
of particles, or the square or the cube of that potential. The topological rules are usual 
Mayer rules, except for some specific rules that avoid double counting (see Section 4.4). 
Screened cluster expansions for other equilibrium quantities are discussed in Section 4.5. 
In particular, formula (|4..31|) allows one to study particle correlations. 

As detailed in Section 5, screened cluster expansions do incorporate physics required 
for treating Saha regime. All Coulomb divergencies are removed via the introduction of 
screened interactions with proper quantum fluctuations. Recombination into complex en- 
tities is accounted for by particle clusters, because truncated Mayer coefficients do preserve 
the non-perturbative structure of genuine Gibbs operators with respect to interactions on 
one hand, while summation over all associated loop configurations restore correct Fermi 
symmetrization on another hand. Moreover, since the screened interaction is close to 
the bare Coulomb potential at short distances, familiar chemical species defined in the 
vacuum do emerge, for instance atoms H, molecules H2, ions H~ and H2 , etc... Statisti- 
cal weights of particle clusters are in correspondance with the densities of the associated 
chemical species formed by recombination of N p protons and N e electrons. Such weights 
incorporate thermal excitations, thermal and pressure dissociation, as well as spectral 
broadening of energy levels by collective effects. Bonds between clusters describe all pos- 
sible interactions between substances present in the gas. 

The present paper is merely devoted to the presentation of the general formalism. In 
Section 5, we also discuss a few possible applications. In forthcoming papers, we shall treat 
quantitatively several physical situations of interest, in particular non ideal contributions 
to Saha equation of state [2S] , the dielectric response of an atomic hydrogen gas [21] and 
van der Waals forces in a partially recombined plasma [2*5] . 
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2 The Hydrogen plasma 



2.1 Definition of the model 

We consider the two-component system (in three dimensions) made of protons and elec- 
trons. In the present non-relativistic limit, the proton and the electron are viewed as 
quantum point particles with respective masses, charges and spins, e p = e and e e = — e, 
m p and m e , a p = 1/2 and a e = 1/2. The kinetic energy operator for each particle of species 
a = p, e with position x takes the Schrodinger form — A where A is the Laplacian 
with respect to x. Two particles separated by a distance r interact via the instantaneous 
Coulomb potential v(r) = 1/r. The corresponding Coulomb Hamiltonian Hmn b for N p 
protons and N e electrons reads 

N n 2 l 

H Np ,N e = + ~'Y j e ai e aj v(\yL i -yL :j \) (2.1) 

i=i ai i^j 

where ./V = N p + N e is the total number of particles. In <j2.1j) . the subscript i is attached 
to protons for i = 1, N p and to electrons for i = N p + 1, ...,N p + N e , so the species index 
on reduces either to p or e while Xj denotes either the position R, of the i-th proton or 
the position rj of the j-th electron (j = i — N p ). 

The system is enclosed in a box with volume A, in contact with a thermostat at 
temperature T and a reservoir of particles that fixes the chemical potentials equal to p p and 
p e for protons and electrons respectively. Its grand-partition function E is (/3 = l/(fceT)) 

E = Tr exp [-0(H Np>Ne - p p N p - fi e N e )] . (2.2) 

In (|2.2j) . the trace is taken over all states symmetrized according to the Fermionic nature 
of each species; the boundary conditions for the wave functions at the surface of the 
box can be chosen of the Dirichlet type. Lieb and Lebowitz [201 have proved that the 
thermodynamic limit (A — > oo at fixed (3 and fj, a ) exists, thanks to Fermi statistics and 
screening. Indeed, the Fermionic statistics of at least one species implies the ^/-stability 

H Np , Ne > -B(N P + N e ), B>0 (2.3) 

that prevents the collapse of the system. On the other hand, screening ensures that it 
does not explode. In a fluid phase, the infinite system maintains local neutrality, i.e. the 
homogeneous local particle densities p p and p e for protons and electrons remain equal for 
any choice of the chemical potentials p a . In other words, the common particle density p = 
p p = p e , as well as all other bulk equilibrium quantities, depend on the sole combination 

p=(p p + p e )/2, (2.4) 

and not on the difference v = (p e — p p )/2. In particular, in terms of the fugacities 
z a = exp(f3p a ), this means that both the density p and the pressure P are functions of (5 
and z = (z p z e ) 1 / 2 = ex.p((3p) only. In the following, we set 

= M ~ ^B^ln^, p e = p + ^rk B T In ^ (2.5) 
A \ p Z \ p 

where X a = ((3h 2 /m Q ) 1//2 is the thermal de Broglie wavelength of species a. This choice 
guarantees that the Maxwell-Boltzmann densities p l p = 2z p /(2ir\ 2 ) 3 / 2 and p l e d = 2z e /(2ir\l) 3 / 2 
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of free (no interactions) proton and electron gases respectively, are identical, i.e. p p = 
pf = 2z/(2vrA 2 ) 3 / 2 with A = (A p A e ) 1/2 . The enforced neutrality 

^e Q z Q /(2vrA 2 ) 3 / 2 = (2.6) 

a 

of the ideal mixture simplifies the analysis of the screened cluster expansions for the 
interacting system derived in Section 4. 



2.2 Formal Mayer expansions 

For describing dilute regimes at low fugacities {z <C 1), Mayer fugacity expansions are 
a priori well appropriate. That expansion for the pressure is easily inferred, at a formal 
level, from identity 

(3P= lj ^, (2-7) 
where it is understood that thermodynamic limit A — > oo is taken once for all. They read 

<"-= £ UK-. p-s> 

where Mayer coefficients B^ p ^ e can be expressed as suitable traces, 

B Np , Ne = j Tr [eM-PH NpiN J] Mayer . (2.9) 
The first Mayer operators [exp(—{3H]y p: ]y e )] Mayer read 

[exp(-/?i/i i0 )]Mayer = exp(-/3H lj0 ) , [exp(-/?# ,l)]Mayer = exp(-/?i? ,l ) , 

[exp(-/3ffi,i)]M«*« = exp(-(3H 1A ) - exp(-(3H 1>0 ) exp(-(3H 0A ), (2.10) 

while similar expressions can be obtained for [exp(— /?i?7v p ,Af e )]Mayer by considering all 
possible partitions of N p protons and N e electrons. The traces (|2.9|) must be taken over 
symmetrized states according to Fermi statistics: for each product of Gibbs operators 
exp(-f3H Mp ,M e ) (M p < N p , M e < N e ) in [exp(-/3#jv p ,jvJ]Mayer, such states are products 
of symmetrized states made with M p protons and M e electrons. For instance, in space of 
positions and spins, -B 2j o reads 



#2,0 = jJdK 1 J dR 2 [4(R 1 R 2 |exp(-/3# 2i o)|RiR2) 



- 4(R 1 |exp(-^ 1) o)|Ri)(R 2 | exp(-/?#i, )|R 2 ) 

- 2(R 2 Ri| exp(-/3i7 2 ,o)|RiR 2 }] (2.11) 

(the prefactors 4 and 2 are due to spin degeneracy). 

Mayer series for particle densities p Pte , are easily obtained by inserting (|2.8|) into the 
identities 

Pa = Zaj^ (2.12) 



for a = p, e. These series read 

CO 

x {n p - iy.N e \ 



00 z n p z N e 

Pv= S (n P --\\\n^ Bn ^ (2 ' 13) 

jV B =lJV e =0 
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and 

00 z Np z Ne 

p*= E N JN e e -i)\ BN *> N °- (2 - 14) 

N p =0,N e =l p y e ' 

Similar Mayer series can be derived for any equilibrium quantity. Contrary to the case 
of systems with short range forces, all these series must be understood in a formal sense: 
Mayer coefficients Bm,n b diverge because the Coulomb potential is not integrable at large 
distances. The first step of our method consists in removing all Coulomb divergencies via 
systematic resummations of chain interactions in quantum Mayer graphs defined for an 
equivalent gas of classical loops (see Section 3). Then, a suitable reorganization of the 
resulting resummed graphs provides the required screened cluster expansion which is well 
suited for describing both low temperature and low fugacity regimes (see Section 4). In 
particular, Mayer series (|2.13|) are rewritten in terms of well-behaved resummed coefficients 
B N N , that do account for recombination into complex entities at low temperature (see 
Section 5). 



3 Path integral formalism 

Within the Feynman-Kac path integral representation, we introduce the equivalent clas- 
sical gas of charged loops (section 3.1). Mayer fugacity expansions for this sytem are con- 
sidered. We briefly sketch the resummation scheme which removes all long-range Coulomb 
divergencies in the Mayer graphs. This provides resummed diagrammatic expansions of 
the equilibrium quantities in terms of a well-behaved screened potential. 



3.1 The gas of charged loops 

In the definition (|2.2j) of the grand-partition function, the trace can be taken over sym- 
metrized Slater sums built with one-body states |x<r 2 } that describe a particle localized 
at x with the z-component of its spin equal to a z . This provides 



00 N p Ne N p N e 



x / Yl dR i II dr i (!) • • (Np)^v e (l) • • - r Pe {Ne) I exp(-(3H Np >Ne ) | R x . . .R^ r 1 ...r Ne ) 

(3.1) 

In (|3.1j) . V Vfi is a permutation of (1, Np >e ), and ^{Vp^) is the signature (±1) of V Pi e- No- 
tice that the spin part of the matrix elements contributes the degeneracy factor * \ frr z \ 

IliJi Ilj=l( a p,Vp(^\ a p,i^K,Ve(J^ a t^ wnich onl y depends on the permutations V p , e . 

The expression (|3.1|) can now be transformed by using the Feynman-Kac path integral 
representation ^3] of all diagonal and off diagonal matrix elements of exp(—[3H]y pj N e ), i.e. 

(x 1 ...x JV |exp(-/?iJ iVp ,iv e )|xi...xjv) = 7- x9 .„,„ — / Vfc) 



exp 



■PHN p ,N e )\*.l-XN) = 11 (27rA2 )3/2 / 11^ 

i=l ^ a i' i=l 

J> ai e Q , f dsv(\(l - s)(xi - x,) + s(x^ - x£) + X^s) - A Q .^(s)|)l (3.2) 

TT. JO J 
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In (|.S.2jl . the functional integration is carried over all Brownian bridges subjected to 
the constraint £j(0) = £j(l) = 0, with the normalized Gaussian measure T>(£) defined by 
its covariance 

' V(^(s)$ u (t) = V inf 0, 0(1 - sup(s, t)) (3.3) 



/ 



In the genuine formulation introduced by Feynman, (1 — s)x- + sx« + \ ai £i(s) is the 
position of the particle i at the dimensionless time s (in units of (3h) along a path joining 
Xj to x^. If we set x- = Rp (j) for i = 1, ...,JV P , and x^ = r-p e (j) with j = i — N p for 
i = iVp + 1, ...,Np + N e , we see that this path is closed if V p (i) = « or V e (j) = j, and 
is open otherwise. However, since any permutation can be decomposed into a product of 
cyclic permutations, all open paths associated with the particles exchanged in a given cycle 
can be collected in a bigger closed path. Thus, as first used by Ginibre [21, 5 then becomes 
identical to the grand-partition function H; oop of a classical system of loops with Maxwell- 
Boltzmann statistics. The Boltzmann factor for loops was reorganized by Cornu ^3] in 
order to exhibit a generalized loop fugacity and two-body interactions between loops. This 
leads to [13 HU 

oo ^ „ N 

S = E loop = E / n^M^n^-W* A)), (3.4) 

^=0 ' i=l i<j 

where the various symbols and quantities are defined as follows. 

First, a loop C is characterized by several degrees of freedom X, a, q, rj(s) with 
s S [0, g]. X is the position of the loop, while a specifies the nature (electrons or protons) 
of the q exchanged particles in the loop. The shape of the loop rj(s) is determined by 
parametrizing the collection of the q trajectories followed by the q exchanged particles as 
X + \aV( s ) with s £ [0, q]. Within this parametrization, the genuine particle positions in 
the matrix elements are of the form x( fc ) = X + \ a rf{k) with k integer, k = 0, ...,q — 1, 
and each trajectory connects x^) to x^ +1 - ) with x^ -* = = X (see Fig.l). The closed 
Brownian path rj(s), tj(0) = rj(q) = 0, is distributed according to the normalized Gaussian 
measure T>(rj) with covariance 

/ B (,M.M*) = Mw(./«,*/«)(i-« ( p(./«,*/«)). (3-5) 

In deriving ()3.5|) . we have used that each individual trajectory is itself a Brownian path 



Fig. 1 - A loop C = (a,q,X,r](s)) 
made of 5 particles. 




X + X a r](s) 



with ending points that are Gaussianly distributed. According to all the external and 
internal degrees of freedom which define the state of a loop, the phase-space measure 
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T>(£) is the product of the discrete summations over a and q, the spatial integration over 
X and the functional integration over rj(s) with the Gaussian measure T>(rj), i.e. 

f v(o... = E E / dX / pfa)- ( 3 - 6 ) 

a=p,e g=l ^ ^ 

The fugacity z{C) reads [T^ l^] 

z(£) = (-r 1 - ^ f 3 exp(-^(/:)), (3.7) 

where the signature (— 1) 9_1 of the cyclic permutation of the q exchanged particles accounts 
for Fermi statistics, the factor 2 arises from spin degeneracy, and the self energy U(C) of 
the loop is 

U(C) = ^ fds f q dt(l-S [s]M )d(s-t)v(\X aV (s)-X a rj(t)\) (3.8) 
z Jo Jo 

with the Dirac comb 

oo 

5(s-t) = E S(s-t-n). (3.9) 

n=— oo 

In (|3.8|) . [s] ([t]) denotes the integer part of s (i), so the factor (1 — #[ a ],[t]) excludes the self 
energies of each individual particle. Eventually, the two-body potential V(£i, Cj) reduces 



V(£ i ,£ j ) = e ai e aj f da [ dt5(s -t)v(\X i + X^s) -Xj - \ a . Vj (t)\). (3.10) 

In (|3. l(Jf) . the Dirac comb H3.9|) ensures that each line element of rj^s) interacts only 
with the corresponding line element of TJjit) taken at a time t which differs from s by 
an integer. This "equal-time" condition can be seen as a manifestation of the quantum 
nature of the particles. Consequently the potential (|3.1U|) between two loops is different 
from the electrostatic interaction energy between uniformally charged wires with the same 
shapes. This difference is the source of algebraic tails in the equilibrium correlations of 
quantum particles [2H1 E0 OOj • 

According to the equivalence formula Q3.4[) , the equilibrium quantities of the quantum 
gas can be inferred from that of the system of loops. For instance, the particle densities 
p a are related to the one-body loops density p(C a ) via the identity [T51 [2*2*] . 

oo „ 

Pa = E / V( Va )q a p(C a ), (3.11) 

where a is the species index of C a . 

The introduction of the gas of loops is quite useful for our purpose because the standard 
methods of classical statistical mechanics can be applied to this system. In particular, the 
equilibrium distribution functions of the loops system are formally represented by series of 
generalized Mayer graphs. In each graph, the usual points are now replaced by loops 

1 Notice that contrary to eq. (5.11) of |22|. the factor e ai e a - is here included in the definition of 
V(Ci,£j). 
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Each loop Ci carries a statistical weight z{Ci). Two loops A and £j are connected by at 
most one Mayer bond 

/(A, A) = exp(-/?V(A,£i)) " 1» ( 3 - 12 ) 

and the other standard topological rules [HI] hold. The contribution of a given graph is 
obtained by integrating over all degrees of freedom of each loop, except some degrees of 
freedom of the root loops that define the arguments of the considered distribution function. 
For instance, loop density p(£ a ) reads, 

P(£a) = ^2^gjz(£ a ) I Y[v(C n )z(C n )[Hf] g . (3.13) 

where the sum runs over all unlabelled topologically different connected graphs Q with one 
root point £ a and N = 0, 1, internal points £ n (for TV" = 0, the integral over the £ n 's is 
replaced by 1). The symmetry factor S(Q) is the number of permutations of the labelled 
loops £ n that leave the product of bonds [T\ f]g unchanged. The first graphs that appear 
in ()3.13j) are shown in Fig. 2. 



° — A A A 

(a) (b) (c) (d) (e) 

Fig. 2 - First Mayer graphs Q . The points stand for loops and the bonds represent 
Mayer factors f = exp(— (3V) — 1. 



Before going further within the loop formalism, it is useful to make the connexion, at 
least at a formal level, between diagrammatic series ([3.13]) and fugacity expansions ([2.13]) . 
Q2.14J1 of particles densities p v ^ e . First, if we sum together contributions of all graphs in 
([3.13]) with a given number N of loops, we recover the Mayer coefficient ,67V that arises in 
the expansion of p(C a ) infered from 

'<« = ^>W- (3 ' 14) 

i.e. 

00 (r \ 

= E (n_1)\ / II T>(£ i )z(£ i )B N (C a ,Ci,... i C N - 1 ). (3.15) 

N=l ^ '' J i=l 

For instance, if the graphs shown in Fig.2a. and 2b. exactly provide B\{C a ) = 1 and 

B 2 (C a ,Ci) = exp(-/3V(C a , C\)) - 1 (3.16) 

respectively, the sum of three graphs with N = 2 shown in Figs. 2c, 2d and 2e, gives the 
third Mayer coefficient 

B 3 {£ a , d, £ 2 ) = exp(-/3V(£ ,£i)) exp(-/3y(£ a , C 2 )) exp(-/3V(£i, £ 2 )) 

- exp(-pV(£ a , A)) - exp(-pV(£ a , A)) - exp(-(3V(£i,£ 2 )) + 2. (3.17) 
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The usual Mayer coefficients -Bjv p ,jv e that arise in expansions (|2.13j) and (|2.14|) of p Pte with 
respect to z p and z e can be recovered by inserting (|3.15|) into the r.h.s. of ()3.11j) . After 
expressing integrations over loops degrees of freedom in terms of that of particles, we 
obtain each coefficient Bn p ,n b by collecting together all contributions with given numbers 
N p and N e of protons and electrons respectively. The expression of Bn p ^ e in terms of the 
suitably symmetrized traces (j2.9|) then follows by applying backwards the Feynman-Kac 
formula 

As mentionned in Section 2.2, all coefficients Bjy p ^ e diverge. Of course, such diver- 
gencies are still present in the Mayer graphs Q. In fact, since the two-body loop interac- 
tions (|3.1()[) are long-ranged and decay as the genuine Coulomb interactions between point 
charges, 

V{£i,£ j )~ ?I ai9j ??. , |X i -X i Hoo, (3.18) 

the spatial integrals over the relative distances |Xj — X,-| between different loops do not 
converge. Such long-range divergencies are easily removed by taking advantage of the 
classical nature of the loops degrees of freedom. The corresponding method is a straight- 
forward extension |161 115j of chain resummations introduced for classical point charges 
|17| I18[ I19| . It is briefly sketched in the next section where we give the resummed form of 
the low-fugacity diagrammatic series H3.13j) for p(C a ). 

3.2 Resummed low-fugacity expansions 

The basic trick which allows us to remove the long-range divergencies consists in the 
resummation of the simple convolution chains built with the long-range part of the Mayer 
bond. The corresponding resummations can be carried out in various ways, according to 
different possible decompositions of the Mayer bond. For our purpose, it is quite convenient 
to make the following choice |20j . 

f(L i ,C j ) = f T (£ i ,L j )-pV(£ i ,£ j ) + ^[V(£ i ,£ j )] 2 , (3.19) 

which defines /r(>Ci, Cj). The replacement of each Mayer bond by its decomposition (|3.19|) 
provides new graphs built with three kinds of bonds, fx, —f3V and f3 2 V 2 /2. Then, all 
simple convolution chains built with bonds — f3V can be resummed in terms of the sole 
screened potential 



4>{Li,Cj) = V(C i ,C j ) + ^2(-l) n (3 n 

n=\ 

/n 
J] V(Ci)z(£i)]V(£ h Cx)...V{C u Ci +1 )...V(C n , Zj) (3.20) 

Once all chain resummations have been performed, the whole set of Mayer graphs Q is 
exactly transformed into a new set of prototype graphs Qp built with the resummed bonds, 

F c = -(3cl) (3.21) 

and 

F R = exp(-/30) - 1 + /30. (3.22) 
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The statistical weights w(C) of loops are either z(C) (bare loop), or z(£)(exp(i#(£)) — 1) 
(dressed loop). Ip is the sum of all simple rings built with —(3V; it is equal to 



I R (C) = \j VifrWdWfaZiWfcuQ = ^(V-<f>)(C,C) (3.23) 

The topological rules which define prototype graphs Qp are identical to those relative 
to Mayer graphs Q, except for the following additional exclusion rules that avoid double 
counting. The convolution of two bonds F c with an intermediate bare loop is forbidden. 
Moreover, when a bare loop is connected to the rest of the graph by a single bond, this 
bond is either F c , or 

FP = F R -?f. (3.24) 

Eventually, the symmetry factors attached to the Gp's are defined as those of the C/'s. 
According to these rules, the whole series (j3.13|) is exactly transformed into 

^)=Ei*) / n^M^n^L' (3 - 25) 

Qp n 

where the sum runs over all unlabelled topologically different prototype graphs Qp while 
w and F denote the various kinds of weights and bonds. The first graphs that appear in 
(|3.25j) are shown in Fig. 3. 



z(J* - 1) 

o owvf o^m c&m 

(a) (b) (c) (d) 

Fig. 3- First prototype graphs Q-p . The possible bonds are F c (diagram b), Fp (diagram 
c) and the ending bond F^ (diagram d). If the weight of a point is not specified, it can 
take the two values z or z(exp(Ip) — 1). 



We stress that the previous remarkably simple result follows from combinatorial and 
topological properties of Mayer graphs (see Appendix A). It does not depend neither on 
the used decomposition of the Mayer bond, nor on the specific form of the two-body poten- 
tial. Therefore, the technical details of these resummations are identical to those already 
displayed in previous works relative to other decompositions of / jl 51 132j or to relativistic 
interactions [SH] (the corresponding proofs are simplified versions of the original derivation 
devised by Meeron ^H] for classical point charges). The decomposition used in ^5] follows 
by expressing the two-body loops potential V as a sum of multipolar interactions. Then, 
the resummations of all convolution chains built with charge-charge interactions lead to 
diagrammatic series analogous to ([3.25ft . The corresponding graphs have the same topo- 
logical structure as the Gp's, with different weights, bonds and exclusion rules. A brief 
comparison between both (equivalent) series can be found in |2f)j . 

The resummed potential <f> is the quantum analogue of Debye potential 20]. In a dilute 
regime, degeneracy effects for free protons and free electrons are weak, so the classical 
Debye screening length kT 1 = (4vr/3 £ Q 4p«)~ 1/2 = (167r/3e 2 z/(2vrA 2 ) 3 / 2 )" 1 / 2 naturally 
emerges. More precisely, at large (but not infinite) distances |Xj — Xj| ~ <p(Ci,Cj) 
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behaves as the classical Debye potential between two point charges q%e. ai and qje a . with 
screening length k -1 , 



At larger distances, |Xj — Xj| S> k^ 1 , l/|Xj — Xj| 3 -dipolar interactions overcome the 
exponentially decaying Debye term. Out of the framework of the loops formalism, this 
mechanism can be seen as an irreducible manifestation of quantum mechanics that ulti- 
mately leads to algebraic tails in the particle correlations. Eventually, at short distances 
|Xj— Xj| <C k^ 1 , <j>(£i,Cj) reduces to the bare potential V(£i,£j) apart from the constant 

According to the previous properties of </>, both spatial and functional integrations in 
any prototype graph Qp do converge j2Hj. The l/|Xj — Xj| 3 -decay of <j) removes all long 
range Coulomb divergencies. Since all Coulomb singularities at the origin are smoothed 
out by functional integrations over shapes of loops, integrability at short distances is also 
fulfilled. This smearing process is nothing but the consequence of Heisenberg's uncer- 
tainty principle which prevents a quantum particle to stay at a fixed position (a similar 
regularization also occurs for Mayer graphs Q built with V). Furthermore, the exponen- 
tiated structure of Boltzmann factors, which is an essential ingredient for the description 
of short-range effects, is preserved through the resummation process. Indeed, both bonds 

(T) 

Fr and provide factors exp(— (3<j)) in the Gp's, which can be replaced by exp(— j3V) 
in a dilute regime: this ensures that recombination can be properly taken into account 
after a suitable reorganization of prototype graphs (see Sections 4 and 5). 

As a conclusion, resummed diagrammatic series 1)3.25(1 is well suited for our purpose. 
The resummed potential <fi does account for classical or quantum mechanisms that prevail 
at either short or large distances. Here, it should be noticed that each prototype graph Qp 
involves infinite sums over numbers qi of particles contained in internal loops Li. Though 
each contribution associated with a given set {qi} is finite because both functional and 
spatial integrations converge, the discrete sum over all q^s may diverge ;2(Jj . Therefore, 
each prototype graph Qp must be understood as a formal infinite collection of graphs, with 
identical structures in terms of loops and different particles numbers inside the loops. 

4 Resummed expansions in terms of particle clusters 

In this Section, series (|3.25j) will be reorganized in terms of new diagrams made with finite 
clusters of particles. In a given cluster, all particles interact together and the weights are 
correctly symmetrized according to Fermi statistics. The reorganization is similar to that 
described in Section 3.1 for recovering the coefficients Bjy p ^ e from the Mayer series (|3.13j) . 

4.1 Introduction of clusters made with fully interacting loops 

In each given graph Qp, if we replace all bonds by their expressions in terms of </>, products 
of Boltzmann factors exp(— (3(f)) appear. A given product of Boltzmann factors then defines 
a cluster of interacting loops. However, such a product may be incomplete in the sense 
that it does not incorporate all Boltzmann factors associated with every pair of loops. 
For instance, if we consider the prototype graph shown in Fig. 4a, the product of bonds 



11^ = F^\C a ,C 1 )F^\C a ,C 2 ) provides the term exp[-/30(£ a , £i)] exp £ 2 )]. 



This term defines the set {£ a , C\, £2} of interacting loops, and it is incomplete since the 
Boltzmann factor exp[— /3<f>(Ci, £2)] associated with pair (£%,£ 2 ) is missing. Since all 




(3.26) 
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Fig. 4 - Prototype graphs showing incomplete 
Boltzmann factors. 





(a) 



(b) 



mutual interactions between particles that belong to loops {C a , C\, £2} are not present, 
the considered graph provides unphysical contributions. Of course, such contributions 
must be cancelled out by other similar contributions arising in other graphs. In fact, if 
we add the contribution from graph shown in Fig. 4b to that from graph in Fig. 4a, we see 
that the previous incomplete product exp[— /30(£ a , C\)\ exp[— j3(j)(C a , £2)] does disappear. 
Here, we first proceed to a systematic reorganization of the whole diagrammatic series 
1)3.25)1 where all incomplete products of Boltzmann factors are removed. 

According to Section 3.2, graphs Gp are identical to Mayer graphs of a system with 
4> in place of V and zexp(In) in place of z, except for some exclusion rules that avoid 
double counting. In the absence of such rules, the (finite) sum of graphs Qp with iV 
loops would provide the Mayer coefficient B^ n, similarly to the derivation of Bn from 
the £7's described in Section 3.1. Since all Boltzmann factors are complete in B^n, the 
required reorganization would then be achieved. Here, the difficulty entirely arises from the 
exclusion rules which prevent to merely perform the previous sum. Our method consists 
in introducing parts of graphs made with M loops which are not affected by the exclusion 
rules, and for which summations over all possible internal structures provide the coefficient 
B(j) : M- The remaining parts are built with either convolutions of bonds F c , or bonds F 2 /2, 
so all Boltzmann factors are eventually incorporated in the Bs m's. 

The first step of the method is based on the introduction of new graphs Qp de- 

(T) 

fined as follows. In each graph Qp, we decompose all ending bonds F R as Fr plus 
—(3 2 (f) 2 /2 = — F c 2 /2. At the same time, we rewrite all dressed weights 2 (£)(exp (/#(£)) — 1) 
of intermediate loops in convolutions F c -kF c as z(£) exp(l#(£)) plus —z{C). The series 
1)3.25)1 then becomes 

^) = E^W / X\ V ^n)w*(C n )\\[F*] (4.1) 

The modified prototype graphs Qp are built with three kinds of bonds F* that are F c , Fr 
or -F c 2 /2, and four kinds of weights w* that are z, z(exp(Ip) — 1), zexp(Iji) or z* = —z. 
They have the same topological structure as the genuine Mayer graphs except for the 
following rules. The intermediate loop in a convolution F c * F c carries a weight which is 
either zexp(l#) or z* . Each bond F 2 /2 is an ending bond that connects a single loop 
with weight z* to the rest of the graph. All other loops (including the root loop C a ) 
carry weights that are either z or z(exp(Ip) — 1). The sum in 1)4.1 J) runs over all non 
topologically equivalent graphs. The symmetry factor S(Q* P ) is defined as usual by virtue 
of the combinatorial properties derived in Appendix A. The first graphs Qp are shown in 
Fig.5. 

In each graph Qp, we can make a partition of the ensemble of loops as follows. On 
one hand, there are loops with weights z* , and on another hand, there remains all other 
loops which can be distributed into clusters. These clusters are defined as the remaining 
connected subsets when the loops with weights z* are removed. In each cluster C, bonds 
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z" 

o OAAt cCm cM&m 

Fig. 5 - First graphs Q v . The ending bond -F c 2 /2 is drawn as a double wavy line. As 
before, the points with no weight specified can take the two values z or zexp(J^) — 1. 



are either F c or Fr, and weights are either z or z(exp(i#) — 1) except for the intermediate 
loops in convolutions F C *F C which carry weights zexp(I^). A given loop in a cluster may 
be connected to another loop, in the same cluster or in another one, by convolution chains 
F c -k F c ... * F c where all intermediate loops carry weights z*; or a similar ring convolution 
with at least two intermediate loops can be attached to this loop. It can be also connected 
to ending loops with weights z* by single bonds F^/2. Several examples of the cluster 
structure of the graphs Q* p are shown in Fig.6. 






Fig. 6 - Cluster structures in the graphs Q v . (a) graphs belonging to the same class 
and leading to a graph Qc when added together, (b) two cluster structures belonging to 
different classes. 



Now we proceed to a classification of graphs Q* p according to their topological structure 
in terms of the previous clusters and connections. All graphs in a given class have the 
same respective numbers of clusters and loops inside each cluster. Moreover, two loops 
connected by convolution chains F c * F c ... * F c (with precribed numbers of intermediate 
loops) belong to the same respective pairs of clusters, or to the same respective clusters; 
similarly, rings F c * F c ... * F c (with precribed numbers of intermediate loops) and ending 
bonds F%/2 are attached to loops which belong to the same respective clusters. Various 
classes of graphs Q* p are shown in Fig.6. In each class, the chains can be attached to all 
the different loops inside a given cluster. Remember that, in such a cluster, the loops can 
carry any weight except z*, and the internal bonds are F c or Fr. Moreover, the connecting 
structure of convolution chains between different clusters are identical. Also, the rings or 
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ending bonds F<?/2 attached to single clusters are identical. 

It turns out that the sum of all graphs Q* p in a given class can be merely rewritten 
as a single graph Qc- Indeed, by virtue of combinatorial properties of Mayer graphs (see 
Appendix A), the various bonds and weights inside each respective clusters, as well as 
the various chains (rings) F c * F c ... ★ F c and ending bonds F^/21, can be added together 
in a simultaneous and independent way. Inside each cluster, the sum of the two possible 
internal bonds F c and Fr provides the Mayer bond = exp(— /3(f)) — 1 associated with <f>, 
while the corresponding sum of weights reduces to z^ = z exp(i^). The further summation 
of all possible internal structures is analogous to that performed when relating usual Mayer 
graphs to Mayer coefficients (see Section 3.1). The corresponding final statistical weight 
W^fC) for a cluster C = L\, Ln made with N loops reads 

wm = Kcec*m Bm (42) 

where BAG) = Bj,,n is the Mayer coefficient for the N loops inside C and interacting via <f> 
(if the root loop L a belongs to C, Nl is replaced by (N — 1)!). The expressions of B^ are 
identical to that of Bn with <p in place of V (for instance see formulas 1)3.16)1 and 1)3.1 7|) 
for N = 2 and N = 3 respectively). The previous simultaneous resummations inside each 
cluster also amount to sum all bonds F c that connect loops belonging to a given cluster C 
to a loop Ci with weight z* . Since F c = — /3(f), this gives rise to the total potential <&(C, Li) 
between C and Li, 

$(C,£) = (4.3) 
while the corresponding bond T c is 

F c = -13®. (4.4) 

Of course, all these summations do not change the topological stucture of the chains, rings 
and ending bonds attached to the various clusters. Eventually, the series 1)4.1)1 is exactly 
transformed into 

= e 5(ky / n »w w *w J n v^z* (L n) [u^] gc [u §} gc ■ 

The integration V{Ci) over the phase space available to cluster Ci reduces to integrations 
over degrees of freedom of internal loops inside Ci, 

V{d) = J] V(L) (4.6) 

The root loop L a necessarily belongs to one cluster, i.e. C a = Co, for which the product 
in ()4.6)) only runs over all the other internal loops in C a . Each bond F c = —(3® connects 
either, one cluster to one loop with weight z* , or two loops with weights z* . All these bonds 
are distributed into convolution chains that connect two clusters, or into rings attached to 
a given cluster. Each bond F^/2\ connects one cluster to one ending loop with weight z*. 
The loops L n with weights z*(L n ) are intermediate loops in the previous chains and rings, 
or ending loops (there are at least two intermediate loops in each ring). Every graph Qc 
is connected. Two clusters may be connected by an arbitrary number of chains, while an 
arbitrary number of rings may be attached to each cluster. The symmetry factor S(Qc) is 
defined as the number of permutations of clusters {Ci,i ^ 0} (except the root cluster C a ) 
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and of loops {C n } with weights z* , that leave the product of bonds [T\ T^g c [Y[ ^c/^] g c 
unchanged. The considered permutations separately act on loops {jC n } with weights z* on 
one hand, and on subsets of clusters with identical numbers of internal loops on another 
hand; permutations of internal loops inside each cluster are already accounted for by the 
1/N\ factor in the definition (|4.2|) of the cluster weight W^(C). Eventually, the sum in 
()4.5|) runs over all topologically different graphs Qc- The first graphs in (|4.5|) are shown 
in Fig. 7. 




Fig. 7 - First graphs Qc, corresponding to the summation of all prototype graphs made of 
one, two and three loops. The second graph in the second row results from the summation 
of the graphs (a) of Fig. 6. The bond T C {C,C) = — /3$(C,£) between a cluster C and a 
loop C of weight z* is drawn as a simple line, while the bond T^jlX is drawn as a double 
line. The statistical weight of a cluster is W^{C). 



In the diagrammatic series IJ4.5J) . the sum of all graphs Qc with a single root cluster 
C a and no attached rings, is nothing but the density of a system of loops with potential 
(f) and fugacity z^: the corresponding expression is identical to (|3.15j) with 4> hi place of 
V and Z0 in place of z. All the other graphs, with at least two clusters or one ring avoid 
double counting: they result from the exclusion rules for prototype graphs Qp generated 
by the chain resummations where V is replaced by <f>. 

As required, each graph Qc no longer involves incomplete products of Boltzmann factors 
exp(— (3(f)). Indeed, all Boltzmann factors exp(— (34>) are reorganized into Mayer coefficients 
13^(0) which define the weights W(p(C) of the clusters. As shown by their construction 
from l|M.14|) . the Bn's involve only complete products of factors exp(— (5V) associated with 
subsets of interacting loops (each product for L loops arises from the L-th order term in 
the expansion of E[ oop with respect to z). Of course, this property is also satisfied by the 
^,jv' s with (j> in place of V. 
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4.2 Introduction of truncated Mayer coefficients 

Similarly to the contribution of any prototype graph Qp obtained by fixing the numbers 
of particles inside loops, the corresponding contributions of every graph Qc are also well 
behaved by virtue of the sufficiently fast decay of the resummed potential <j>. In a dilute 
regime, at finite distances, (f> tends to V, so B^n can be merely replaced by Bn in weights 
W^(C). However, this replacement is no longer legitimate for large distances between the 
loops inside cluster C: it is crucial to keep screening effects embodied in eft for ensuring 
integrability. In order to circumvent this difficulty, we introduce truncated Mayer coeffi- 
cients B^ N such that their analogues Bjj with V in place of 4> remain integrable at large 
relative distances. Of course, there are several ways of defining truncated coefficients by 
appropriate substraction of terms which are not integrable when (j) is replaced by V. The 
definition used below is aimed at preserving the topological structure of Mayer graphs. In- 
deed, within this definition, we find that the resulting diagrammatic series has a standard 
Mayer form in terms of the Bjj's (see (|4.18|0 , 

Of course, no truncation is necessary for Ba,,\ = 1 since there is only one loop in the 
associated cluster. So, we first consider B^p that is given by formula ()3.16j) with <p in place 
of V. We define B^ 2 as 

^,2 = ^,2 + /^--^- + -^-- (4.7) 

By construction, B^ 2 behaves as /3 4 (/> 4 /4! at large distances. Therefore, B 2 (£i, £2) decays 
as l/|Xi — X2I 4 when |Xi — X2I — * 00, so B 2 indeed is integrable over the whole space. 
Notice that integrability at zero separation results from the smearing process already 
evocated in Section 4.2. 

In order to define B^ 3 , we replace each Boltzmann factor exp(— f3(j)) into B^^ by 

ex P (-/30) = 1-/30 + + Bl 2 . (4.8) 

The resulting expression for B^ s(Ci, £2, £3) is a (finite) sum of terms that is nothing but 
the sum of Mayer diagrams where the Mayer bond has been decomposed according to 
(|4.8|) . This provides various bonds g between two loops which can be either F c , F?/2\, 
F c 3 /3! or BT 2 . Then, B^[C\, £2, ^3) can be rewritten as 



B^(£ 1 ,£2,£3) = Y J [U9} D3 - ( 4 - 9 ) 



In l|4.9|) . L>3 is a connected labelled diagram made with the three loops (£1,^2,^3), such 
that two loops are connected by at most one bond g. Moreover, [n#]D 3 denotes the 
product of bonds in D%, and the sum is taken over all the different diagrams (for which 
the products [rjjjflj are different); of course no integration is performed in D3 over 
loops degrees of freedom. The truncated coefficient B^ 3 is defined as the contribution to 
(|4.9jl of all diagrams which remain integrable with respect to all relative distances when 
cj) is replaced by V. For diagrams which are built with bonds F c , F c 2 /2! or F c 3 /3!, the 
integrability condition is fulfilled by a simple power counting. Indeed, since V decays as 
l/|Xj — Xj|, the sixth-dimensional integration over the two relative distances converges if 
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the overall power of bonds F c is larger than, or equal to 7. This gives 

+ E + E B IA?% + E ^,2^ 2 f 

7712 772 7712 7713 7713 773 

+ Z^^.2 2! 2 , +Z^^2 2! 3 , +Z^ £ V>,2 3 , 3, 

^ c 3! 3! ^ 2! 2! 3! ^ 2! 3! 3! 3! 3! 3! ' 1 ' ' 

where the sums run over all the (labelled) diagrams D3 with the considered products of 
bonds £>J 2 , F c , i ? c 2 /2! and Fj}/3\. According to this definition, we rewrite (|4.9|) as BT 3 
plus the remaining terms that are not integrable when <j) is replaced by V: 

B w = ^,3 + E B h F c + E B h F ^ + E ^ff + E B l^% 

773 772 773 772 772 

p2 773 773 773 772 173 

+ E tt + H 3f t + FcFcFc + S FcFc ir + E ^ t 

F 2 F 3 p2 p2 p2 p2 p2 

+ E Fc iJif + E Fc ^if + • ( 4 - n ) 

The sums have the same meaning as in (|4.1()|) in terms of diagrams D%. 

Now we introduce new diagrams Dq by identifying clusters made with one or two loops 
in the diagrams D%: two loops are in the same cluster if they are connected by a bond 
jBJ 2 . In these diagrams Dc, the loops of D3 are replaced by the corresponding clusters. 
The weights of clusters with one and two loops are B^ ; % = 1 and 23 J 2 respectively. Two 
clusters are linked by at most one of the three possible bonds T c , f^/2\ or T^/3\, where 
J- c is related to the total interaction potential <£ between the clusters, i.e. T c = — /33>. The 
total potential $ between a cluster with two loops and a cluster with one loop is given 
by (|4.3|) . while it obviously reduces to <j) if two clusters are made with single loops. The 
bonds between a two- loops cluster and a one- loop cluster, result from summing together 
the contributions of all diagrams D3 where the two given clusters are linked by bonds F c 
for T c , by bonds F C F C or F c 2 /2! for J*/2\, by bonds F c F c 2 /2! or F c 3 /3! for ^ c 3 /3!. Similarly 
to (|4.11|) . the resulting expression for B^ can be detailed as 

b« = s« + E + E B lM + E B l^r 

■r-2 773 773 773 772 773 

+ E^f : 3 f + Ef ^+^^c+E^^+E^if 

772 773 772 772 772 772 772 

+E^ifif+E^if# + ififif- <" 2) 

where each sum runs over all labelled diagrams £>e (loops inside a cluster are labelled) 
with the considered product of weights and bonds. We stress that definition 1)4.10(1 allows 
to introduce clusters which are connected by single bonds T c , T^/2\ or T^/3l. If B^ 3 
had been defined as the sum of diagrams with bonds £>T 2 only (the first two terms of 
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PTDI V left-over terms in (1X171 like E ^2^,2^ or E ^2 F c /2!^ c 2 / 2! > could not have 
been rewritten in terms of clusters connected by the sole links J- c , J^j1\ or 

The construction of the general truncated coefficient B^ N follows from a procedure 
similar to that introduced above for B^ 3 . Each Boltzmann factor is replaced by (|4,8|) into 
B<t>,N, with the result 

^=E[ruL (4 - i3) 

n N 

where diagrams Djy are defined with the same rules as the D^s. We define B^ N as 
the contribution to (|4.13|) of all diagrams -D/v that remain integrable with respect to all 
relative distances when (f) is replaced by V. For diagrams built with bonds F c , or 
F c 3 /3!, the simple power-counting rule used for £^3 implies that the overall sum of powers 
of bonds F c is strictly larger than (3iV — 3). However, for N > 4, this rule is necessary 
but not sufficient. For instance, the diagram shown in Fig. 8 fulfills the counting rule, but 
it is not integrable (when <fi is replaced by V) with respect to all separations. In fact, 
the determination of the integrable diagrams requires a closer analysis of their topological 
structure. 



Fig. 8 - A diagram that fulfills the power counting rule but 
that is not integrable with respect to the relative distance [R4 — 
R3I when <j) is replaced by V. This diagram does therefore not 
appear in the definition of B^ 4 . It rather belongs to the class of 
graphs characterized by a F^-link between the two clusters of loops 
{£i,£ 2 ,£ 3 } and {£4}. 



In each of the non- integrable diagrams, we define clusters made with M loops (1 < M < 
N), for which integrals over relative distances between loops inside a given cluster do 
converge when <f> is replaced by V. By construction, these clusters realize a partition of 
the ensemble of N loops of the considered diagram Djy. As a result of this procedure, 
two connected clusters are linked by either one bond F c (i^-link), two bonds F c or one 
bond F c 2 /2! (F c 2 -link), three bonds F c , or one bond F c and one bond F^/2\, or one bond 
F^/3l (i^-link) . A diagram Dq is defined as the sum of all diagrams Dn with identical 
partitions into clusters and identical structures of F™-links (1 < n < 3) between clusters. 
For a given cluster made with M loops, the summation of all possible internal structures 
of bonds g provides, by definition, the truncated coefficient B^ M (<6</,,i = 1 if M = 1). 
The bond T$ between two connected clusters, which results from the sum of all possible 
F c n -links, reduces to T c for n = 1, T 2 c j2\ for n = 2, or T?/3\ for n = 3. Now, T c = 
where $ is the total interaction potential between two clusters, i.e. 

m,Cj)= E E <p(£kXi),- (4.14) 

The detail of the summation leading to T^j2\ is illustrated in Fig. 9. 
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Fig. 9 - Diagrams D4 belonging to the same class characterized by a F^-link between 
the two clusters of loops {£i,£2} &nd {£3, £4} (the bond B^ 2 is drawn as a gray link). 
Their sum leads to the V c diagram BT 2 (£i, £ 2 )BT 2 (£ 3 , £ A )^/2\ in (l4~T51) . 



Therefore, similarly to formula (|4.12|) for N = 3, we can express Ban i n terms of truncated 
coefficients as 

b*n = ®In + e [ n b Im] d [ n ?*] D ■ ( 4 - is ) 

D c C C 

In Q4.15JI , the sum runs over all labelled and connected diagrams Dq with at least two clus- 
ters (all M's are strictly smaller than N), while [II^Afl-Dc anc ^ H^l-Dc denote, respec- 
tively, the product of weights and the product of bonds for a given diagram Dc- Moreover, 
this sum is restricted to the diagrams which become not integrable over all relative dis- 
tances between clusters when cp is replaced by V. For instance, in the explicit expression 
(|4.12j) of (j4.15|) for N = 3, the diagrams Dq with three one-loop clusters and products of 
bonds F c F*/3\F3/3\, ^/2\J^/2\^/3\, J* f2\ J* or J^ c 3 /3!^/3! J*/3\, do not 
appear. When ./V increases, a detailed account of the forbidden diagrams Dc becomes 
rapidly cumbersome, like the explicit evaluation of the truncated coefficient B^ N itself. 
However, we stress that this difficulty does not affect further diagrammatic reorganizations 
which are only based on the general structure of (|4.15j) . 

Now, we return to sum (J4.5|) over graphs Qc- In each graph, we insert, into the 
definition (|4.2|) of each statistical weight IV^(C), the expression (|4.15|) of Mayer coefficient 
B(f,(C) = B^^n f° r a cluster C made with TV loops. The whole series 1)4.5(1 is then rewritten 
in terms of new graphs which are still made with clusters of loops. The weight Wj(C) 
of a cluster C is 

w t {c) = UceeM^) B T {c) (4 16) 

Like in (|4.2|) . for root clusters C a which contain the root loop £ a , N\ is replaced by (N — 1)1 
in (|4.16|) . Two clusters are connected by at most one bond T which reduces to either T c , 
J^/2\ or For any one-loop cluster C = £ which is doubly connected to the rest of 

the graph by two bonds JF C , weight (|4.16|) is replaced by 

Wj(£) = z(£)(eMW)) ~ !)• (4-17) 

This specific weight results from the summation of two contributions W^C) = z^(£) = 
z(£) exp(Iji(£)) and z*{£) = —z(£) associated respectively with a convolution T c * T c in 
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a diagram Dc, and from the same convolution in a graph Qc- Similarly, for a one-loop 
cluster C = C which is singly connected to the rest of the graph by one bond .F c 2 /2!, the 
weight WT(C) also reduces to (|4.17|) . Therefore, series (|4.5[) is exactly transformed into 




(4.18) 



where the sum runs over all unlabelled topologically different graphs Qq. Like (|4.15|) . 
this sum is restricted to graphs which become not integrable with respect to the relative 
distances between clusters when <fi is replaced by V. The symmetry factor S(Qq) is the 
number of permutations of labelled clusters that leave the product of bonds [n^>]gg' 
unchanged (only clusters with identical numbers of loops are permuted, except the root 
cluster Co = C a ). Like in (|4.5|) . the permutations of internal loops inside a given cluster 
are accounted for by the factor 1/iV! in definition Q4.16JI of cluster weight Wj(C). This 
rule follows from the same combinatorial properties (see Appendix A) that have been 
used in the resummation process of section 3.2 or in the construction of graphs Qc- Here, 
such properties allow us to transform the sums (|4.15|) over labelled clusters (with labelled 
internal loops) into sums over unlabelled objects. The first graphs in (|4.18f) are shown in 
Fig. 10. 




Fig. 10 - A few graphs Qq with one or two clusters. The bonds T$ between two clusters 
can be either T c = — /?$ (single line), y\T1 (double line) or Jj^c (triple line). The small 
while point represents root loop C a . The clusters are drawn as shaded circles comprising 
a certain number of loops. The statistical weight of a cluster is given by the truncated 
weight WT(C) KWh . 



As required, each cluster weight WT(C) remains integrable with respect to relative 
distances between its internal loops when <j) is replaced by V. Notice that products of 
Boltzmann factors embedded in Wj(C) are complete: because of the structure of decom- 
position (|4.15j) . the truncated coefficients N , like the B^jv's, do not involve incomplete 
products. Moreover, the introduction of specific weight (|4.17f) can be traced back to the 
exclusion rules for prototype graphs Qp that avoid double counting. 
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4.3 Introduction of clusters made with particles 



If we insert (|4.18|) into the expression (|3.11j) of proton density p p , we find 

-. oo „ „ n 

pp=T,sm^ / p (^)^/n^x^)[n^] e -' (4 - l9) 

gT V C ) q a= \ i=o C 

where root loop C a inside root cluster C a in each graph Qg, is made with q a protons. When 
the phase space integrations over loops inside each cluster C% are explicited, according to 
()3.6|) . in T>{Cj) given by (|4.6|) . each graph Qq provides an infinite sum of contributions 
where clusters C(M p , M e ), made with M p protons and M e electrons, appear. In general, 
in each individual contribution, the weight associated with particle cluster C(M p , M e ) 
is not correctly symmetrized according to Fermi statistics. For instance, the graph in 
Fig. 10c. provides, for (a a = p,q a = l;a± = e,q% = l;a% = e,q2 = 1), a cluster C(l,2) 
with a weight that does not include any exchange between the two electrons. However, 
the graph in Fig.lOb also provides, for (a a = p,q a = l\a\ = e,q± = 2), the same cluster 
C(l, 2): when both contributions are summed together, Fermi symmetrization is restored. 
In the following, we rearrange the whole series in terms of new graphs G made with particle 
clusters that incorporate Fermi statistics. 

Let us consider the class of graphs Qq made with the same set of loop clusters {C{} 
and the same topological connecting structure of bonds Two graphs in a given class 
differ by the number of internal loops in at least one cluster Cj. For instance, the graphs 
in Figs. 10b. and 10c. belong to the same class defined by a single root cluster C a . On 
another hand, the graphs shown in Figs.lOd., lOe. and lOf. belong to different classes 
because respective bonds T§ between clusters C a and C\ are different. For each cluster C, 
specification of the nature and of the number of particles in all internal loops, leads to the 
introduction of a particle cluster made with N p protons and N e electrons. In each given 
class, once all graphs Qq have been replaced by infinite sums over loop particle numbers, 
we sum together the contributions associated with identical sets of particle clusters {Cj} 
(all respective numbers (A^ , iV- ) are identical), connected by the same bonds J-a. This 
defines a graph G. 

In a given graph G, each particle cluster C(N p , N e ) can be traced back to a loop cluster 
C with ./V loops. In the sums ^ obtained by expliciting V(C), there are several terms 
which can be associated with particle cluster C(N p , N e ). Each of them involves L p protonic 
loops (a.i = p) and L e = N — L p electronic loops (ctj = e), while the corresponding ordered 
sets of particles numbers, [q[°\ <Tl \ with q[ a ^ > q^ > ... > q^\ are such that 

Y^ql a) =N a . (4.20) 
k=l 

The internal state of a particle cluster C then appears to be defined by the partitions 
Q a = [gj , qj* ] of N a , and by the internal degrees of freedom of L a loops with particles 
numbers determined by Q a , L a being in the range 1, ...,N a . The corresponding cluster of 
loops is denoted by C(Q p ,Q e ). An example of all possible partitions Q a and clusters of 
loops C(Q p ,Q e ) associated with the same particle cluster C(3, 2), are shown in Fig. 11. 

The statistical weight Zj(C) depends on internal state (Q p ,Q e ',C(Q p ,Q e )) of C. In 
the sum £\ arising from T>(C), there are N\/{L p \L e \) terms which define a loop cluster 
C(Q P , Q e ) made with L p protonic loops and L e electronic loops. In the corresponding sum 
^ , there are L p \L e l/(Y\^L 1 n p (q)ln e (q)\) terms associated with the same partitions Q p 
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Particle cluster C(3, 2) 


Qp Qe 
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[1,1,1] 


[1,1] 





Fig. 11 - The six possible partitions (Q p , Q e ), and the associated set of loops C(Q P , Q e ), 
involved in the internal state of the particle cluster C(3, 2) made of 3 protons and 2 
electrons. 



and Q e , where n a {q) is the number of loops with q particles of species a. Therefore there 
are 

JL X W = N} (421) 

different terms in Yl ai wn i cn contribute to Zj(C). Once integrations over positions 
and shapes of loops are performed, all these terms provide identical contributions. Thus, 
we define ZT(C) as the product of the number of terms (|4.21|) by the weight Wj(C(Q p , Q e )) 
of loop cluster C(Q p ,Q e ). The phase space measure T>(C) involves all possible internal 
states of C(N p , N e ): it reduces to a discrete sum over all partitions (Q p , Q e ) satisfying the 
total particle number constraints (|4.2()j) . combined to spatial and functional integrations 
over positions and shapes of loops inside C(Q P , Q e ). 

The common connecting structure of graphs Qq remains unchanged through the sum- 
mations leading to graphs G. Two particle clusters Cj and Cj are connected by at most 
one bond T$ which can be either —/?<&, /3 2 $ 2 /2!, — /3 3 <3? 3 /3!; here $ is the total interaction 

potential between two loop clusters Cj(Q^, Q\ ) and Cj(Q ( j'\ Q^) that describe internal 
states of Ci and Cj respectively. Again, by virtue of the combinatorial properties derived 
in Appendix A, S{Q^) is exactly transformed into S(G) which is the number of permu- 
tations of clusters Cj (except C a ) with identical particle numbers (N^ p \ N- e ^) that leave 
the product of bonds unchanged. Similarly to Q4.19JI . the resulting series only involve 
graphs G which are no longer integrable over the relative distances between the particle 
clusters when cj) is replaced by V. 



4.4 The screened cluster expansion 

According to section 4.3, diagrammatic series (|4.19|) is exactly transformed into the so 
called screened cluster expansion, i.e. 

^ = E^) / V(C a )ZT(C a )q a J \{V{C l )Zl(C l )[\{^] G (4.22) 
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The graphs G are identical to the usual Mayer graphs, where the points are now particle 
clusters, except for some specific rules. The rules which define the graphs G are listed as 
follows: 

• Every graph G involves the root cluster Co = C a and n (n > 0) black clusters (Ci, ...,C n ). 
The root cluster C a contains the root proton which is fixed at the origin in p p . Each cluster 

Cj (i = 0, ...,n) contains N^' protons and electrons. 

• The internal state of a cluster C(N p ,N e ) (C G {Ci,i = 0, ...,n}) is determined, on one 
hand, by partitions Q a = [q[ a \ qj_ a ] with q[ a ^ > > ... > q£ (a = p, e); on another 

hand, by loop cluster C(Q P , Q e ) made with L a loops ) where carries q^ 

particles of species a (a = p,e). The numbers n a (q) of loops carrying q particles of species 
a satisfy the total particle number constraints 

N a 

Y / qn a (q)=N a . (4.23) 

q=l 



• The phase space measure T>(C) for a cluster C(N p ,N e ) reads 

V (°) = E / ft dX i P) ft dX l £) / ft ft V(r,i e) ) (4.24) 

Q P ,Qe k=l k=l k=l k=l 

where xj^ and denote position and shape of loop respectively. The discrete 
sum is carried out over all partitions (Q p ,Q e ) satisfying ()4.23j) . For root cluster C a , loop 

df^ is identified to root loop C a {q^ 1 = q a ), and no integration is performed over position 

(p) 

= X a which is fixed at the origin. 

• The statistical weig ht Zj(C) for a cluster C(N p ,N e ) reads 

= U h? iC fl Z *ft ) EF + m » Qe)) - (425) 

1 1<j=1 n p\Q)- 11<j=1 n e{q)- 

For root cluster C a , n p (q a )\ is replaced by (n p (q a ) — 1)!. For a cluster C different from 

(p) 

C a , the internal state of which is determined by a single protonic loop D[ or a single 
electronic loop c[ , the expression (|4.25|) is replaced by (a = p, e) 

ZT(C) = z,(c[ a) ) - z(4 a) ) (4.26) 

when C is either, the intermediate cluster of a convolution J-^J-a,, or connected to the rest 
of the graph by a single bond J r ^ ) /2\. The truncated Mayer coefficient £>T N is defined by a 
suitable truncation (see Section 4.2) of the usual Mayer coefficient jS^jv f° r N loops with 
pair interactions 4>. This truncation ensures that BT n remains integrable over relative 
distances between loops when <p is replaced by V. The first truncated Mayer coefficients 
are 

fl£i = l, ^ 2 = ex P (-W-l + ^-^ + ^f, ... (4.27) 
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• Two clusters Cj and Cj are connected by at most one bond J~^{C^ Cj) which can 
be either -/?$, /? 2 $ 2 /2!, -/3 3 $ 3 /3!. The potential $(Ci,Cj) is nothing but the total 

interaction potential between loop clusters Ci{Q\ , q\ ) and Cj(Qj,Qj) that describe 
internal states of Cj and Cj respectively, i.e. 

SKC^Cj) = $(C i5 C,) = E ^(A^)- (4-28) 

The product of all bonds inside G is denoted by H The graph G is connected, 

i.e. for any pair (Cj, Cj) there exists at least one connecting path made with one or more 
bonds Tfy. 

• The symmetry factor S(G) is the number of permutations of labelled black clusters Cj 
(i > 1) that leave the product [n^^G unchanged. Only clusters with identical numbers 
of protons and electrons are permuted. 

• The sum ^ G runs over all topologically different unlabelled graphs G which are no longer 
integrable over relative distances between clusters {Ci,i = 0, ...,n} when <f> is replaced 
by V. 




Fig. 12 — A few graphs G occuring in the screened cluster expansion of p p . (a) graph made 
of a single proton; (b ) graphs involving a proton and an electron; (c) graph consisting in the 
cluster C(3, 2), whose possible internal states are listed in Fig. (d) further examples. 



Several examples of graphs G are shown in Fig. 12. We stress that graphs G have the 
same topological structure as familiar Mayer diagrams: the ordinary points are replaced 
by particle clusters, while usual Mayer links are now bonds Ta,. The introduction of 
specific weight ()4.26|) avoids double counting arising from elimination of V in favor of 
(j). For a system with short range forces, a diagrammatic series similar to ()4.22|) can 
be derived without any prior chain resummations: (j) ls merely replaced by the genuine 
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N p ,N e 



potential, and all the above rules are conserved except that associated with specific weight 
(|4.26[) . Of course, in this case, the sum of all graphs with the same number N of particles 
can be rewritten in terms of the sole standard Mayer coefficient of order N. Here, the 
(finite) sum of all graphs Gm p ,n e with the same total numbers N p of protons and N e of 

electrons, provides a contribution, that can be rewritten as (zp p z^/{N p - l)\N e \)B {R) 

by definition of renormalized Mayer coefficient N ^ . Then, screened cluster expansion 
(|4.22j) can be recast as 

oo Np TV 

*= E w^^U- <«»> 

Np = l,N e =0 y p ' 

The long range divergencies which plague the bare coefficient B^ p ,N e , (see Section 2.2) are 
removed: every graph Gn p ,n e is finite because $(Cj,Cj) decays as a dipolar interaction 
for large separations between Cj and Cj. Since screened potential <f> depends on fugacities, 

( R} 

B N N does not depend on the sole temperature, contrary to usual Mayer coefficients for 
systems with short range forces. 

4.5 Particle correlations and other observables 

The particle correlations, as well as other physical observables, can be inferred from multi- 
body loop densities. Hence, the screened cluster expansion of such physical quantities will 
directly follow from the expansions of these loop densities that are analogous to (|4.22|) . In 
particular, by starting from (|4.18|) and following the same lines as in Section 4.3, we easily 
obtain for (one-body) loop density p{C a ) 



I G 



G ^ ' J CeCa J i=l 

(4.30) 

where graphs G are identical to those in (|4.22|) . The sole difference between (j4.3U|) and 
(|4.22j) lies in the integration upon degrees of freedom of root cluster C a . Now T>(C a ) and 
Zj(C a ) are also given by (|4.24j) and Q4.25J) respectively, while factor 5(C,£ a ) selects the 
contribution of internal states of C a where one loop C is identified to C a . It is easily 
checked that series (|4.22j) for the protonic density can be recovered by inserting (|4,3()|) 
into ()3.11|) . Notice that factor l/(n p (q a ) — 1)! occuring in 1)4.22(1 results from the fact that 
there are n p (q a ) loops in (|4.30|) identified to C a . The screened cluster expansion of other 
one-body physical quantities might be also merely derived from (|4.3(J|) . 

The truncated two-body loop density (loop correlations), can be expanded similarly 
to (|4,3U|) . Its screened cluster expansion reads 



PT{CaM = J2s(G) [ V{Cab) ( £ 5{C k ,C a )5{C h C b ))zl{C ab ) 
G { } J c k ,dec ab ,k^i 

i 

n v / J r, r ]eC b 

Y[V(Cd2$(C i )[l[r+\ (4.31) 



x 
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Now, graphs G involve, either a single root cluster C a b that contains both loops C a and 
Cb, or two roots clusters C a and C& with C a in C a and in C&. The rules that define such 

graphs with n black clusters {Cj} are identical to that of (|4.22jl . In particular, the graphs 

(p) 

are connected, and only permutations of black clusters Cj with identical numbers iV- 

and are considered in the definition of symmetry factors. In 1)4.31)1 . the measure and 
weights associated with root clusters are identical to that of black clusters, as in 1)4. 3U|) . 



5 Physical interpretation of graphs G and concluding re- 
marks 

Expansion 1)4.29)1 is a reorganization of formal series (|2.13[) that is suitable for studying 
Coulomb gas at low fugacities and low temperatures, whether in a fully or partially ionized 
phase. Indeed, graphs in that expansion admit a natural interpretation in terms of the 
recombined entities and their mutual interactions, as discussed below. 

• Graphs with a single cluster 

First, we consider graph Gn p ,n e with a single root cluster C a . When fugacities go to zero, 
(ft reduces to V, and we can replace all ring factors exp(/fj) by 1 and all Mayer coefficients 
B^ by B T , in Zj(C a ). Then, integration over all internal states of C a reduces, after 
applying backwards Feynman-Kac formula, to the trace of a suitably truncated Mayer 
operator [exp(— /^A^A^Mayer- F° r shown in Fig.l3a, using expression (|4.7|) of <Bj 2 , 
we find 

[expt-ZHTi,!)]^ = eM~PH ltl ) - exp(-/?^ 1 ) ) 

+ / driexp[-(/3-ri)i?f 1 ) ]yi,iexp[-riiI 1 (0 1 ) ] 
Jo 

- I dn r dr 2 exp[-(/3 - t^H^V^ exp[-(r 1 - t 2 )h[°}]V 1>1 exp[-r 2 H[ }] 
Jo Jo 

+ tin dr 2 / dr3exp[-(/3-r 1 ) J ffJ L ) ]n,iexp[-(r 1 -r 2 )ilJ 1 ) ] 
Jo Jo Jo 

x Vi,i exp[-(r 2 - T 3 )H[°I]V 1A exp^tfg], (5.1) 

where h[°1 = H\ q + H$\ is the kinetic part of while Vx t % is its potential part. 

Similar expressions for [exp(— /3i/7v p ,7V e )]Mayer can De obtained from the definition of *Bj N 
(see Section 4.2). In the resulting truncated structure analogous to 1)5 each Gibbs factor 
is associated with a complete Hamiltonian Hm p ,m £ {M p < N p and M e < N e ) that does 
incorporate any pair interactions between M p protons and M e electrons. At low fugacities, 
the leading contribution of Gn p ,n b to B N Ng then reads 

Truncation in [exp(— /3HN p ,N e )]wia.yer> inherited from that in the B T, s, ensures that trace 
1)5.2)1 is finite contrary to that of [exp(— /^ifjVp.ArJlMayer defining B^ p: ]y e : when it is ex- 
pressed in position and spin spaces, matrix elements of [exp(— /3#7Vj,,7V e )]Mayer are inte- 
grable for any spatial configuration, i.e. their decay is sufficiently fast with respect to any 
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large particle separation. Moreover, that trace is correctly symmetrized with respect to 
Fermi statistics as a direct consequence of summation in T>(C a ) over all possible 

distribution of particles into loops. 

According to the above properties, trace ()5.2|) does account for recombination into complex 
entities. Indeed, when temperature goes to zero, Gibbs operators associated with Hm p ,m e 
provide exponentially growing terms arising from possible bound states with negative 
ground state energies E$ M . Single term (M p = N p , M e = N e ) can be interpreted as the 
ideal contribution, of the possible complex entity made of N p protons and N e electrons, 
while product of terms (M p < N p , M e < N e ) describe its dissociation into several entities. 
Ideal contributions associated with familiar entities can be singled out in various graphs. 
For instance, G\ : \ shown in Fig. 13a provides the ideal contribution 

Az p z e [(mp + m e ) /{2tt (3h 2 )] 3/2 exp(-pE H ), (5.3) 
of atoms H in their groundstate, or Cr2,2 shown in Fig. 13b provides the ideal contribution 

2z 2 p zl [(2m p + 2m e )/(2irph 2 )] 3/2 exp(-/3E H2 ) (5.4) 

of molecules H2 in their groundstate with energy E H2 = Ef\. Notice that trace fl5~2"]) 
appears as the natural (N p + iV e )-extension of the so-called second virial coefficient first 
introduced by Ebeling [2] for N p = N e = 1 (in that latter case, leading behaviour (|5.3j) of 
Q5.2j) can be readily inferred from perturbative low-temperature expansions |3|). 
If leading contributions to the present graph Gn p ,n b arise from simple or complex entities 
in their ground state with bare interactions, next corrections involve thermal excitations 
and collective effects. For instance, in G2,2> corrections to (|5.4j) arise from excited states 
of molecule H2, as well as from thermal dissociation into atoms H, ions H2 and H~ , 
protons and electrons. First corrections due to screening are obtained by expanding, in 
ZT(C a ), all ring factors exp(i^) in powers of Ir and all Mayer coefficients B^ in powers 
of f3(4> — V) (both Ir and (3(4> — V) go to zero in a sufficiently dilute regime). At finite 
fugacities, non-perturbative collective effects like pressure dissociation, are accounted for 
in £>J. Since <j> decays faster than the Coulomb potential, excited bound states with typical 
sizes larger than the screening length disappear: for instance, Rydberg states of atom H 
with principal quantum number n, are removed from contribution of G\ t \ (see Fig. 13a) for 
n 2 as > k^ 1 . Moreover, spectral broadening of energy levels in vacuum, directly follows 
from double-time dependence of (j)(d,£j) in Sj and Sj (any line element of £j interacts 
with any line element of Cj). Since the equal-time condition that defines V (see (|3.10|) ) 
no longer holds for <ft |2l)| . functional integration of cannot be expressed in terms of an 
effective Hamiltonian. This is a particular manifestation of a general property regarding 
the motion of quantum particles coupled to a thermalized bath |34| . 

• Graphs with two clusters and more 

According to the previous analysis, a particle cluster C(N p ,N e ) can be associated with a 
possible complex entity made with N p protons and iV e electrons. Therefore, a graph G 
with several clusters describes interactions between such entities. We consider successively 
various contributions arising from bonds T${Ci,Cj) between either neutral or charged 
particle clusters C\ and Cj. 

Let us consider two neutral clusters Cj and Cj (numbers of protons and electrons in each 
given cluster are identical). Even if screening is incorporated into ^(Cj, Cj), at sufficiently 
low fugacities ^(Q, Cj) can be replaced by the bare interaction between neutral clusters. 
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Indeed, like eft itself, that interaction displays a dipolar structure at large distances, and 
its contributions to G remain integrable. For instance, in graph Cr2,2 shown in Fig. 13c, 
collective effects in Z$(C a ), Zj(Ci) and *(C«,Ci) can be neglected in a sufficiently di- 
lute regime: the corresponding leading contribution accounts for familiar van der Waals 
interactions between two atoms H. 

Now we turn to a neutral cluster C, connected to a charged cluster Cj by a bond /3 2< 3? 2 /2! 
or — /3 3 $ 3 /3!. At sufficiently low fugacities, $(Cj, Cj) can be replaced, by the bare dipole- 
charge interaction between Ci and Cj, the square or cube of which is indeed integrable 
at large distances. For instance in graph £2,3 shown in Fig. 13d, leading contribution 
describes polarization effects of an atom H submitted to the electric field created by an 
ion H~ . In a bond — /3#, screening is a priori crucial for ensuring integrability (however, 
for some graphs sufficiently connected, leading contributions can still be estimated by 
keeping only bare interactions). 

Eventually, we consider pairs (Cj, Cj) of charged clusters with a connecting bond T${Ci, Cj). 
Effective interaction <&(d,Cj) in each bond takes into account screening due to free pro- 
tons and free electrons at a mean-field level. Suitable combinations of charged pairs and 
bonds in G, provide contributions that account for various mechanisms. For instance, 
graph G^i shown in Fig.l3e gives screened contributions beyond mean-field Debye-like 
terms. Graph shown in Fig.l3f describes modification of screened interactions be- 
tween free protons due to the presence of ions H% ■ Summation of all similar graphs with 
convolution chains (—/?$) * (—/?$) * ... * (—/?#) and intermediate charged complex entities, 
allows us to introduce effective interactions screened by all charges present in the medium. 
Polarization of charged complex entities by other charges also appear. Graph G^i shown 
in Fig.l3g incorporates polarization effects of an ion by a free proton which is itself 
screened over distances of order k^ 1 . 

Thus, screened cluster expansion (|4.22j) provides a complete framework that sustains 
the so-called chemical picture at sufficiently low temperatures and low densities. Any 
complex entity made with N p protons and N e electrons can be associated with a particle 
cluster C(N p , N e ), while both ideal and interaction contributions are incorporated into 
statistical weights Zj(C) and bonds Every physical mechanism (see examples above) 
is described by a well behaved graph G, and consequently it can be estimated in a perfectly 
controlled and unambiguous way. We stress that screened cluster expansions ([4.22J) , (|4.30[) , 
Q4.31J1 can be extended to any quantum nucleo-electronic plasma. Then particle clusters are 
made with any kind of nuclei and electrons, while statistical weights and bonds account for 
specific masses, charges and spins of all nuclei species. All topological and combinatorial 
rules that define graphs G, remain unchanged. Usefulness of screened cluster formalism for 
studying partially ionized gases is illustrated by various applications as described below. 

It has been rigorously shown |35U36| that the equation of state of the Hydrogen plasma 
reduces to that of an ideal mixture of protons, electrons, and atoms H in their ground 
state, in some suitable zero-temperature scaling limit defining the so-called Saha regime. 
The first rigorous result on the atomic limit has been established by Fefferman [3Zj (for a 
review, see jZZl). Screened cluster expansions allow us to recover immediately ideal terms 
associated with free protons, free electrons and atoms H. All other contributions are shown 
to decay exponentially faster in the Saha regime. For instance, the ideal contribution (|5.4I) 
of molecules H2 is smaller than (|5.3|) by factor exp[{3(2fi+ Eh — Ejj 2 )\ ~ exp[/3(3i£# — Ejj 2 )] 
(for /i close to Eh), which indeed exponentially vanishes at zero temperature since 3Eh < 
Eh 2 - though molecule H2 is more stable than two atoms H (Eh 2 < 2Eh), entropy 
dissociation overcomes energy minimization in the present very dilute regime. Moreover, 
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Fig. 13 - Examples of graphs G describing various physical effects. Formation of hound 
states: (a) hydrogen atom (b) hydrogen molecule H2; (c) van der Waals interactions 
between two H atoms; (d) polarization of atom H by the field of a H~ ion; (e) contribution 
beyond mean-field to the screened interactions between two free (ionized) protons and a 
free electron; (f) modification of the screened interaction between two free protons due 
to a H2 ion; (g) polarization of an ion fL^ by the held of a free proton screened over 
distances of order k^ 1 . 



detailed analysis of graphs G provide exact estimations of various corrections to ideal 
leading terms at finite (low) temperature and density. Such corrections arise from excited 
bound states of atoms H, thermal and entropy dissociation, recombination into other 
complex entities (H2, i?^, H~ ,...), screening, interactions between entities, diffraction 
and exchange (all entities behave as almost classical particles). Those results will be 
published elsewhere [23] . 

For the Hydrogen plasma in Saha regime, particle correlations have been determined 
[25] by using screened cluster expansion 1)4.31(1 . That analysis leads to the introduction 
of effective interactions between any entity, which all decay as 1/r 6 at large relative dis- 
tances r. In particular, familiar van der Waals forces between atoms H do appear. At zero 
temperature, standard perturbation formula for two atoms in the vacuum is exactly recov- 
ered. Corrections to that formula due to excited states and screening, are also calculated 
at finite temperature and density. 

The dielectric response of the Hydrogen plasma in its atomic phase has been inves- 
tigated with screened cluster formalism |38) I24| . That phase is defined within previous 
Fefferman limit |37| where \x is chosen slightly larger than Eh'- it reduces to an ideal 
mixture of atoms H in their groundstate (such phase may also be viewed as a particular 
Saha regime where the fraction of ionized protons and electrons is arbitrarily small). The 
behaviour of wavenumber-dependent dielectric constant is analyzed by using expansion 
()4.31j) in linear response to an infinitesimal external charge. A cross-over phenomenon 



31 



emerges, from a dielectric behaviour at wavenumbers much larger than k, to a conductive 
behaviour at scales smaller than k. Conductive behaviour is due to a finite tiny amount 
of free charges that are always present at finite temperature and density (here, screening 
length k _1 is much larger than the mean distance between atoms H). The large plateau 
observed at scales larger than k, leads to a natural definition of the (finite) dielectric con- 
stant for the atomic phase, independently of the presence of residual free charges. At zero 
temperature, the usual expression in terms of the polarizability of a single atom H is then 
exactly recovered. 

Previous examples confirm that screened cluster formalism is well suited for studying 
dilute regimes at low temperature. This concerns both conceptual problems for Coulomb 
matter, like the existence of van der Waals forces at finite temperature and finite density or 
the definition of a dielectric regime, and also derivation of exact asymptotic expressions for 
equilibrium quantities. All those various analysis can be applied to any partially ionized 
(even neutral) gas. In addition, screened cluster expansions provide a promising route 
for deriving either approximation schemes or effective models. Both schemes and models 
might reasonably describe regimes with higher temperature and density, where analytic 
estimations are no longer possible. 

A Combinatorial properties of Mayer graphs 

We derive various combinatorial properties for Mayer graphs which are used in the trans- 
formations described in Section 4. These properties are related to the following central 
observation. Let G be an unlabelled Mayer graph built with N internal points x, weights 
z and bonds / (for our purpose it is unnecessary to take explicitly into account the root 
points). Its contribution reads 



where S(G) is its symmetry factor, i.e. the number of permutations of the ./V internal 
points that leave unchanged the product of bonds [Q fij\Gi ah associated with a labelled 
diagram Gi a b with the same topological structure as that of G. The number of different 
labelled diagrams Gi a b is equal to N\/S(G). Once all integrations over internal points 
(x±, xn) have been performed, all contributions from different labelled diagrams are 
identical. Thus, ()A.1|) can be rewritten as 



where the sum runs over all different labelled diagrams G[ a b with the given topological 
structure determined by G. The transformation of (|A.1|) into (jA.2|) allows us to derive 
several useful identities as described below. 

• Decomposition of bonds 

If bond / is decomposed into a sum, / = g + h, contribution (jA.lj) can be rewritten 
in terms of that of graphs G* built with two kinds of bonds, g and h. Replacing each 
bond fij by gij + hij into (|A.2|) . we find a sum over labelled diagrams G* lab where the 
product of bonds reads [n^fc/lG^ b times [Q h mn ]c* ■ Obviously labelled diagrams G*[ ab do 
exhibit the topological properties of Mayer graphs: they are connected, and two points are 







32 



connected by at most one bond which is either g or h. All GJVJs have different products 
ill Stole?? jn^mn]<?,* b - They can be distributed into different classes where each class 
is defined by an unlabelled graph G* (a = a,b, ...). In each class, two different labelled 
diagrams Gf b can be inferred from each other by a permutation of labels of internal points, 
so there are N\/S(G^) such different labelled diagrams G*[ ab which all provide the same 
contribution after integration over the cc's. Thus, we obtain 

^/[n d -(-)] G [n/]c=E^/[n d -(-)]cjn^c S [n /i ] G , ( a - 3 ) 

where the sum runs over all different unlabelled Mayer graphs G* a built with two kinds of 
bonds g and h. Moreover S"(G* ) is the usual symmetry factor defined as the number of 
permutations of internal points that leave the product [I]to]G* b [Y\^mn}G* b unchanged 
for any labelled diagram G* lab associated with G* . The topological structures of graphs 
G* are inherited from that of G: they only differ by number and position of bonds g and 
h. If we sum (|A.3jl over all possible Mayer graphs G, we then find the remarkable identity 



G* 



(A.4) 

where both graphs G and G* are built within the same Mayer prescriptions. 

• Decomposition of weights 

Similarly to (|A.4f) . we can transform ()A.1|) when weight z is decomposed as z = u + v. 
Again in (|A.2|) we replace each z; L by Ui + Vi. Following the same lines as above, we obtain 
the identity 

E^y/tn <m*)] a tn /] a = e / n d *] n g* n «&] g* n /] 

(A.5) 

where now Mayer graphs G* are defined with two kinds of weights u and v (and one kind 
of bond /). The symmetry factor S(G*) is defined as the number of permutations of 
internal points that leave products of weights [n^ilc??,,! [Il^ilc??, and bonds Uito]©?, 
separately unchanged. Notice that an identity analogous to (|A.5|) holds when the internal 
state of x is made of two components s and t, i.e. dx = ds + dt: u(x)dx and v{x)dx are 
merely replaced by z(s)ds and z(t)dt respectively in the r.h.s. of l|A.5|) . 

• Resummations of internal structures 

Now, let us consider labelled diagrams Gi a b associated with the following Mayer graph G. 
Weights can be either u or v, while bonds can be either g, h or h? /2!. Two connected 
points with weights u are linked by a bond g. A point with weight v is either singly 
connected to a point with weight u via a bond h 2 /2l, or is the intermediate point of a 
convolution h*h (the extremal points of convolutions h * h * ... * h carry weights u). The 
formula analogous to that relating (|A.1|) to (|A.2|) reads here 

^/[n d 4[n^] G [n^] G [i]4[n^[n^] G = 

N 

wT, n d ^[n^c !a jn^ Gia jn^] Gia jnM Gia jn^/ 2! ] Giab - ( a - 6 ) 



1=1 
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Then, we define a partition of points x with weights u in terms of M clusters X, such 
that the subdiagram associated with X is connected via bonds g only. The set of clusters 
X is determined by the topological structure of the unlabelled graph G. Each cluster Xi 
contains Ni points x with weights u, so there remains L = N — ^2i=\ N% points x with 
weights v. There are iV!/(iVi!iV2!...iVA/!L!M!) different partitions of N labelled points into 
M labelled clusters (Xi,X%, ...,Xm) and L labelled points with weights v. Such labelled 
clusters and points (weights v) together with bonds g and h, define a labelled diagram 
Qiab- The corresponding N\/ '(JV1IiV2L.iVfcfI.LlM!) labelled diagrams Gi a b provide identical 
contributions in the r.h.s of (IA.6I). Their sum reduces to 



■ J i=i i ' 1 k=i 

(A.7) 

np\Giab * s a product of convolution chains h * h * ... * h, with intermediate 
points Xk (weight v{xk)) and two ending points belonging either to the same cluster Xi or 
to two different clusters Xi and Xj. Moreover, in [f\ hg r /2\]g lab , each bond h 2 /2\ connects 
a point Xk (weight v(xk)) to a point belonging to a cluster Xi. Now, we can proceed to a 
further summation of contributions (|A.7|) which only differ by the points (inside clusters 
Xi) at the extremities of convolution chains h * h * ... * h or bonds h 2 /2!. This obviously 
leads to the introduction of 

H(X u x k )=^2h(xf,x k ), (A.8) 
i=i 

so (|A.6|) can be rewritten as 

isjfE/n [ n T^n^],,n^)[n^] a jn^] Sl . l - 

Slab »=1 fc=l 

(A.9) 



In l|A.9|) . [Y\ H np ]g lab is the product of convolution chains H * H * ... * H (intermediate 
points Xk and H(xf.,xi) = h(xk,xi)) that either are attached to one cluster Xi or connect 
two clusters Xi and Xj. Moreover, in H 2 r /2\]g lab , each bond H 2 /2l connects a point 
Xk to a cluster Xi. The sum runs over all labelled diagrams Qi a b with different products of 
weights and bonds, and the same topological structure defined by G. Identity (|A.9j) can 
be easily summed over all graphs G that determine the following set of labelled diagrams 
Qiab'- clusters Xi, points xp., bonds H np and H 2 r /2l, are identical, but internal structures 
of clusters in terms of bonds g are different. Since combinatorial prefactors for each 
labelled diagram Qi^ are all identical, the sum over labelled diagrams of products (over 
i = 1,...,M) of internal weights [Jlffimlxj reduces to the product (over i = 1,...,M) of 
sums of those weights over all possible internal structures. This provides for each cluster 
Xi the weight Zi 

where the sum runs over all different products of bonds gi m connecting labelled points 

(i) 

inside X{, i.e. over all corresponding Mayer labelled diagrams built with N points x) 



34 



and bonds gi m . If we set 

Ni 

dX i = Y[dx < ?\ (A.ll) 
j'=i 

we eventually find 

_ i /• M L 

=t,um\] n^™n^^*)[n^]ft Ot [ii fl Sr/ 2, ] « 

Glai) »=1 fc=l 

=E^y/ [n^wi^n^^tn^gtn^!], (a.i 2 ) 

where the last line follows by applying backwards the relation between unlabelled graphs 
and labelled diagrams. The permutations accounted for in S{Q) act separately on sets 
of labelled clusters (X\, ...Xm) and labelled points (x\, xl), while they keep invariant 
each specific product of weights or bonds (this implies that permuted clusters necessarily 
contain the same number of particles). 

All present identities result from simple combinatorics rules on one hand, and from the 
topological invariance of Mayer graphs under the considered transformations on another 
hand. They apply for any kind of points x: usual points (without internal degrees of 
freedom), loops (extended objects), clusters (sets of entities). Results derived in Section 
4 (and those recalled in Section 3.2) follow by successive applications or combinations of 
these identities. 
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